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I. INTRODUCTION 



The expansion of our universe has been one of the greatest attractions of scientific talents since the seminal work 
of Edwin Powell Hubble 1] in 1929. Hubble's compilation of observational distance-redshift (expressed in terms of 
radial velocity) data suggested a linear pattern of "extra-Galactic nebulae" (an archaic term for galaxies) receding 
from each other: 



where H is the proportional constant now bearing his name, and x is the positional coordinates of a galaxy measured 
with our Galaxy as the origin. 

The discovery of Hubble's Law marked the commencement of the era of quantitative cosmology in which theories 
of the universe can be subjected to observational test. Since the days of Hubble, advances in technology have enabled 
astronomers to measure the light from increasingly deeper space and more ancient time, and our ideas of the entire 
history of the expanding universe have been gradually converging into a unified picture of Big Bang-Cold Dark Matter 
universe. In this picture, the dominating form of energy density transited from radiation to dark matter, and relics 
of primordial perturbation were imprinted on today's observable CMB anisotropy and large-scale structures (LSS). 
This picture is obtained from its two ends: the CMB last-scattering surface at z « 1000 and the LSS around us at 
z rs 0. The vast spacetime extent between both ends, in particular the era before reionization, remains mostly hidden 
from our view. In addition, the past two decades' cosmological observations, especially those of type la supernovae 
(SNIa), indicated that the recent history of universal expansion is an acceleration, possibly driven by an unknown 
"dark energy" [2j, |3| whose physical nature has not been identified. 

Therefore it appears to us that our understanding of the universe is currently under the shade of three dark clouds 
— the mysterious dark energy that drives late-time accelerated expansion, the nature of dark matter that is vital to 
the formation of structures, and the unfathomable dark age that has not yet revealed itself to observations. This is 
the "3-D universe" in which possible answers to some of the most profound questions of physics are hidden. 

In the face of these vast unknown sectors of the universe, any observational probe into its past history is invaluable. 
Recently, the direct measurement of the expansion rate, expressed in terms of the Hubble parameter H(z), is gaining 
increasing attention. As a cosmological test, it can help with the determination of important parameters that affects 
the evolution of the universe, and reconstruct the history around key events such as the turning point from deceleration 
to acceleration. As an observable, it manifests itself in various forms in different eras, especially in the baryon acoustic 
oscillation (BAO) features in the LSS that may be detectable in the dark age. 

This paper is a review on the current status of observational Hubble parameter data and its application in cosmology. 
In Section [II] we briefly review the cosmological background of an expanding universe. In Section IlIII we present two 
important observational methods of H(z) observation, their principles and implementations. Next, we review the 
important role of the observational H(z) data in the study of cosmological models in Section HVl We will also discuss 
some issues associated with their application. Finally, in Section [Vj we briefly discuss some ongoing efforts that 
promise possible improvements over the current status of H(z) measurements. 



In this section, we will review some basic ideas and definitions in cosmology that must be kept in mind in order to 
understand and interpret the observational H{z) data and their implications. 



The spacetime structure of the homogeneous, isotropic, and spatially flat universe is characterized by the Friedmann- 
Robertson- Walker (FRW) metric 



x = Hx, 



(1) 



II. BACKGROUND 



A. Spacetime, Metric, and Coordinates 




(2) 



The presence of the scale factor a(t) means that the spacetime is not necessarily static. In reality, we know that the 
universe is expanding, and a(t) increases with time. 
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Using the metric ([2]), the infinitesimal spacetime interval scalar ds 2 = dx^dx" = g^dx^dx" is obviously 

ds 2 = -c 2 dt 2 + a 2 (t)oVoV. (3) 

Here we have used the four-coordinate vector x a — (ct, x 1 ) that has the dimension of length. 

It is often useful to express the spatial components of the four-coordinate vector, i.e. the "comoving position" , in 
dimensionless spherical coordinates x % = (r, 8, <j>) in order to extend the metric to non-flat situations, and give the 
scale factor the dimension of length. Under this convention, the spacetime interval ([3|) can be re-written as 

ds 2 = -c 2 dt 2 + a 2 (t) ( + r 2 d6 2 + r 2 sin 2 ddcj) 2 ) (4) 
\ 1 — kr J 

where k is one of {—1, 0, 1}. The parameter k is the sign of the spatial curvature, and k = if the universe is spatially 
flat. 

We can further transform equation (|4]) by introducing the coordinate 

r dr ' ■ -i m 

Jo v 1 — kr' z 
where the sinn function is a shorthand notation: 

fsina; for k = 1, 
x for k = 0, (6) 

sinha; for k = — 1. 

Switching to the spatial coordinates (x, 8, 0), the interval ds 2 can be written as 

ds 2 = -c 2 di 2 + a 2 {t) [d X 2 + sinn 2 X (d0 2 + sin 2 0d(j) 2 )] . (7) 

The physical interpretation of \ can bee seen by placing ourselves at the origin r — and consider a distant, 
comoving photon emitter in our line-of-sight direction with the coordinate r = r e . Rotate the coordinates so that the 
direction of the emitter has 8 = 0, 4> = 0, we find 

d, 2 Wd £ 2 + a 2 (t )T i^ (8) 

along the line-of-sight. Let t e be the time of photon emission and to that of its reception. Since light-like worldlines 
have ds 2 = 0, we find, for the photon: 

*° cdt r* dr . . 

X(r e ). (9) 



a(t) Jo Vl - kr 2 

Consider the integrand in the left-hand side of equation (J9)). The line element da; = cdt is the physical distance the 
photon has traveled during the time interval dt. But by dividing the physical distance by a(t) we get the comoving 
distance, therefore \ can be interpreted as the total, integrated comoving distance between the emitter and us. If the 
space is flat, this comoving distance is just the difference in the radial coordinate Ar — r e — = r e . 

Sometimes it is convenient to introduce the conformal time, or the comoving horizon rj as the time component of 
the four-coordinate. The conformal time is defined as 

^-Lw) (io) 

where we integrate from the "beginning of time" . Using cq as the time component, the comoving four-coordinate can 
be written as a dimensionless vector x a = (cry, \, &, </>) and the FRW metric takes the form 



9p» = a 2 {n) 



/-I \ 
1 

sinn 2 x 

V sinn 2 x sin 2 8 J 



(11) 
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B. Expansion, Redshift, and the Hubble parameter 



In the introduction we mentioned Hubble's Law discovered in 1929. Hubble's original paper had profound impact 
upon the history of astrophysics and, to a greater extent, mankind's perception of the universe, but here we only take 
some time to appreciate two of his timeless insights. 

At the end of his paper Hubble briefly discussed the possible mechanisms for "displacements of the spectra" 
(i.e. redshift, in modern terms) in the de Sitter cosmology model in which the expansion of the universe is dominated 
by a vacuum energy. He pointed out the two sources of the redshift: the first being "an apparent slowing down of 
atomic vibrations" and the other attributed to "a general tendency of material particles to scatter" . In today's words, 
the first is the special-relativistic effect of Doppler shift caused by the peculiar motion of galaxies, and the latter the 
general-relativistic, cosmological redshift which is linked to the expansion of the comoving grid itself. In the rest of 
this article we will see how these two effects arise in modern cosmology and end up in our observational figures. 

Hubble also noted that his proportional law might be "a first approximation representing a restricted range in 
distance", therefore deviating from the pure de Sitter model in which the Hubble constant H should indeed be 
constant everywhere and throughout the history. This is exactly how we see it now. In the contemporary context, we 
usually define the Hubble parameter H to be the relative expansion rate of the universe: 

H=~, (12) 
a 

and its value is usually expressed in the unit of km s _1 Mpc -1 . The Hubble constant, H , now officially refers to the 
current value of the Hubble parameter. 

However, it is not apparent how this definition is related to observable quantities. Therefore we have to relate 
equation (TT2")) to physical observables such as the length, the time, and the redshift. 

First, we note that the cosmological redshift z at any time t is related to the scale factor a. Let t e be the time of a 
photon's emission by a distant source and to the time of its reception by an observer "here and now." The observed 
redshift z of the source satisfies 

1 + *=^. (13) 

a{t e ) 

Consider an observer who surveys various sources with different redshifts. The ideal survey is assumed to complete 
instantly — all the observations are done at exactly the same time instance to. Of course this is not strictly true, 
but we do not expect the scale factor a(to) to change "too fast", and we expect the redshift not to change too much 
during the temporal scale of our interest (i.e. typical lifetime of humans or observation programs). If we do allow 
t to change however, we are led to the Sandage-Loeb test [3, [H] that observes the drifting of redshift during a long 
period of time. Recently, the variation in the apparent magnitude of stable sources over to has also been proposed 
as a possible cosmological test @. To our best knowledge, no data have been produced using these methods by now, 
and the proposed observation plans usually require ~10 years to yield meaningful results (however, we note that 
the idea of "real-time cosmology" is gaining interest recently, as reviewed by Quercellini et al. @). In this paper we 
will not focus on these methods, and we therefore neglect the passing of to. 

We therefore differentiate equation (fl~3|) with respect to t e , setting to as a constant: 



da(t e ) _ a(t ) dz _ a(t e ) dz 
dt e ~ (l + z) 2 dtT ^^1 + zdtT' ^ ' 

Dividing both sides by a(t e ) we immediately find 

, n 1 dz , _. 

^ = -TTIdV (15) 

In Section [ill Al we will see how equation (|15[) is useful in measuring H{z) by observing passively evolving galaxies. 

Another way to relate H(z) to observable quantities is to use the notion of the comoving distance x introduced in 
equation ([5]). Take the time derivative of equation ©, we find 

dt e a(t e ) 



On the other hand, equation (|14[) tells us about another derivative dt e /dz. Therefore we can find the derivative of \ 
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with respect to the redshift: 



dx = dx dt e 
dz dt e dz 



c a(t e ) dt e 

a(t e ) (1 + z) da(t e ) 

ca(t e ) dt e 

a (to) da(t e ) 
c 



a(t )H' 



that is, 



d [a(to)x\ 
dz 



H(z) 



(17) 



(18) 



(also see, for example [10|,lll|, but beware of different notation conventions). If an observable object spans the length 
a(t )Ax along the line-of-sight in some redshift slice Az, we can estimate H{z). But how do we find such objects, 
i.e. "standard rods"? The idea is not to use the length of a concrete object. Instead, we explore the spatial distribution 
of matter in the universe and focus on its statistical features, such as the BAO peaks in the two-point correlation 
function of the density field. This is another method for extracting H(z) data from observations. (The quantity 
a(to)x can be seen as a distance measure. It is closely related to the "structure distance" ds = a(to)r defined by 
Weinberg [l2|, Chapter 8] that naturally arises in calculating the power spectrum of LSS. From equation ((SJ we can 
see that the structure distance is equivalent to a(to)x if the space is flat, or if the object is not too far away.) 

We remark that the derivation of H(z) expressed in terms of the standard rod, equation (|18p. is only part of the 
story, for we have only considered a standard rod placed in the line-of-sight direction. The transversely aligned test 
body is related to another important cosmological measure, namely the angular diameter distance Da{z) — a(z)r(z). 
In an expanding universe, the angle A9 subtended by a distant source is 



Ad = 



a{z) 



At i = 



a (t ) 



iX + z)D A {z) 



Ar±, 



(19) 



where Ar± is the transverse spatial span of the source measured in the difference of comoving coordinate r [H, [l4| ■ 
Naturally, once the physical scale of BAO is known and the BAO signal measured, the corresponding angular diameter 
distance can also be used as a cosmological test. 

A classical cosmological test is the Alcock-Paczyhski (AP) test [l5| that can be expressed as another combination 
of H (z) and Da(z). The observable of the AP test is the quantity A(z) — Az/(zA8) of some extended, spherically 
symmetric sources, where Az is the difference in redshft between the near and far ends of the object, and A8 the 
angular diameter. By our equations (fT8|) and (fl9|) it can be expressed as 



A(z) = 



Az 

1~a1 



l±±D A (z)H{z)gL. 
z Ar± 



(20) 



A well-localized object placed in a region not too far away from us (so the non-trivial spatial geometry can be 
neglected) will have A% ss Ar n , the difference in the comoving coordinate along the line-of-sight. Furthermore, for a 
nearly spherical object the approximation Ar n ~ Arj_ holds, and A(z) is reduced to 



A(z) = 



l + z 



D A {z)H{z) 



(21) 



Clearly it cannot constrain H(z) or D A (z) separately, but a combination of both. The AP test, in more modern 
context, is usually understood as a geometrical effect on the statistical distribution of objects instead of concrete 
celestial bodies (see [l6l - [l8| . and also OH HI] where the BAO effects were explicitly treated in the analysis). 

Another combination of H (z) and Da (z) naturally arises in the app lication of BAO scales measured in the spherically 
averaged galaxy distribution, namely the distance measure Dy [20J defined by 



D v (z) 



cz(l + z) 2 D 2 A (z) 
H{z) 



1/3 



To break the degeneracy between H(z) and D A {z) in Dy(z), the full 2-dimensional galaxy distribution must be used, 
with the correlation function conveniently decomposed into the line-of-sight and transverse components (see section 
IIIIB1 but also see [2l| for another decomposition scheme). 
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III. HUBBLE PARAMETER FROM OBSERVATIONS 

Equations (fT5j) and (p~8|) are the bare-bone descriptions of two established methods for H(z) determination: the 
differential age method and the radial BAO size method respectively. Either has been made possibly only by virtue 
of state-of-the-art redshift surveys such as the Sloan Digital Sky Survey (SDSS) 22]. In this section, we will review 
both methods and the data they produced. 

A. The Differential Age Method 

As equation (|15l) suggests, to apply age-dating to the expansion history, we look for the variation of ages, At, in 
a redshift bin Az [23|. The aging of stars serves as an observable indicator of the aging of the universe, because 
the evolution of stars is a well-studied subject, and stars' spectra can be taken and analysed to reveal information 
about their ages. However, at cosmological distance scales it is not practical to observe the stars one by one: we 
can only take the spectra of galaxies that are ensembles of stars, possibly of different populations. Since different 
star populations are formed at drastically different epochs, it is important for us to identify galaxies that comprises 
relatively uniform star populations, and to look for more realistic models of star formation. 

The identification of such "clock" galaxies and the observation of their spectra have been carried out for archival 
data [H, and surveys such as the Gemini Deep Deep Survey (GDDS) H, VIMOS-VLT Deep Survey (VVDS) and 
the SDSS [26} . In addition, high-quality spectroscopic data have been acquired from the Keck I telescope for red 
galaxies in galaxy clusters [27j |. Among the galaxies being observed, special notices should be paid to the luminous 
red galaxies (LRGs). LRGs are massive galaxies whose constituent star populations are fairly homogeneous. They 
make up a fair proportion in the SDSS sample and, beyond serving as "clocks" , also trace the underlying distribution 
of matter in the universe (albeit with bias). Therefore, they reveal BAO signature in the density autocorrelation 
function that is used as the "standard rod" in the size method. 

The identification and spectroscopic observations of these galaxies have led to direct determinations of H (z) in low 
and intermediate redshift ranges. Jimenez et al. [24[ first obtained a determination of H{z) = 69 ± 12 km s^ 1 Mpc -1 
at an effective redshift z w 0.09 by the differential age method. The work was later expanded by Simon et al. [28j 
who extended the determination of H(z) to 8 more redshift bins up to z ~ 1.8. This dataset was brought up-to-date 
by Stern et al. [H, Table 2], Recently, new age- redshift datasets for different galaxy velocity dispersion groups have 
been made available (2i| from SDSS data release (DR) 7 LRG samples. We will see how these data are used in the 
study of cosmology models in Section ITVl 

One may wonder why we take the effort to calculate the age differences in redshift bins when the age (or lookback 
time) data themselves can also be used to test cosmological models. Indeed, the absolute age has been very useful in 
the estimation of cosmological parameters (30l - [33| . Nevertheless, precise age-dating with low systematic biases can be 
only carried out on a narrow selection of sources. On the other hand, by taking the difference of the ages in narrow 
redshift bins, the systematic bias in the absolute ages can hopefully cancel each other [34|. Of course, we are not 
gaining anything for nothing even if the systematics perfectly cancel, for the binning of data lowers the total amount 
of measurements we can have. 

A further approximation is that the majority of stars in the galaxies are formed almost instantaneously, in a single 
"burst" [35| , therefore the intrinsic spread of the measured age arising from a heterogeneous star formation history can 
be expected to be small when fitting the observed spectra to stellar population models (specifically the single-stellar 
population (SSP) model used in [28] and [26]). However, recent developments in the study of the formation history of 
galaxies and their stellar populations have led us to re-consider the assumptions made in previous works. For example, 
using galaxy samples selected from numerical simulations, Crawford et al. [36| have shown that the SSP assumption 
may contribute to the systematic bias that varies across redshift ranges (hence failing to cancel, and propagating into 
the differential ages), while models that take the extended star formation history into account can be used to reduce 
the errors on H(z). 

In addition to the complexities in the stellar populations in each galaxy, the heterogeneity of galaxies in the sample 
also contributes to the errors in H(z) measurements. In [361 ] . new sample selection criteria have been proposed that 
could help with obtaining more homogeneous galaxy samples for future analyses. 

B. The Radial BAO Size Method 

In Section III B[ we mentioned that the "standard rod" we seek in the sky is not an actual object but a statistical 
feature. Indeed, the physical sizes of distant celestial objects are usually poorly known. Worse still, even the apparent, 
i.e. angular, sizes of galaxies are ambiguous because galaxies do not show sharp edges, and they appear fuzzy in images. 
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It can be imagined that size measurements along the line-of-sight could only lead to more problems, because even the 
angular sizes cannot help us much in this case. Therefore, identifying a statistical "standard rod" becomes a necessity. 

In the study of LSS, correlation functions are a simple and convenient measure of the statistical features in the 
spatial distribution of matter in the universe. (For an early yet important treatment of the topic in the context 
of galaxy surveys, see [37j- For an example of other statistics in the context of BAO, see |38j.) The two-point 
autocorrelation (i.e. the correlation of a density field with itself) function £(ri,r2) is one of the most used member 
in the correlation function family. It measures the relatedness of position pairs in the same density field: the joint 
probability of finding two galaxies in volume elements dVi and dV2 located in the neighborhood of spatial positions 
j"i and T2 respectively is 



where n is the mean number density. If we believe that our universe is homogeneous in a statistical sense (i.e. that the 
probabilistic distribution, or ensemble, from which the densities anywhere in our particular instance of the universe 
is drawn, does not vary from one area of the universe to another), the autocorrelation function becomes a function 
of r = r*i — r-i only. If we further assumes the (statistical) isotropy of the universe, the direction of r becomes 
unimportant, and the autocorrelation is dependent on the magnitude of r only (that is, £ = £(?"))■ Actually, our 
assumption of homogeneity is unnecessarily strong if we only work with two-point statistics, and all we need is the 
homogeneity in the first two moments of the underlying ensemble. Such an ensemble is known as a wide-sense 
stationary (WSS) one. 

For a WSS ensemble, the famous Wiener-Khinchin theorem says that the autocorrelation and the power spectrum 
P(k) form a Fourier transform pair: 



(Here we write the power spectrum as P{k), independent of the direction of the wave vector k, under the same 
assumption of statistical isotropy mentioned above, but see discussion about redshift distortion below.) Therefore, 
either the power spectrum or the autocorrelation can serve as a statistical tool to reveal the information contained 
in the LSS. Methods of estimating P(k) has been developed and the importance of the power spectrum emphasized 
[39l Eoj . On the other hand, for BAO surveys the autocorrelation function is probably a more straightforward way 
of presenting the results and testing their significance, because the BAO scales manifest themselves as protruding 
features ("peaks" or "bulges") in £(r). Actually, an estimator to the autocorrelation, along with its variance, can also 
be conveniently constructed from survey data using pair counts between the survey and random fields [4l| . 




Needless to say, the "true" autocorrelation of the ensemble can never be fully known, because we have only one 
realization of the random field which is the universe we live in. However, estimating the autocorrelation still makes 
sense because for today's large and well-sampled surveys the assumption of ergodicity is valid, under which the 
statistics can be performed to infer knowledges about the underlying ensemble |12l Chapter 8 and Appendix D] . 

Thus, if a random process induces some features in the spatial distribution of matter, the autocorrelation can be 
numerically computed to reveal such features that are otherwise hidden in the seemingly stochastic distribution. Fur- 
thermore, if the mechanism and properties of this process is well understood and quantitatively modelled, parameter 
estimation using these features becomes a possibility. 

One of such possibility is provided by the BAO signatures in the LSS. The mechanism of BAO effects must be 
traced back to the early universe before recombination, when the Compton scattering rate was much higher than the 
cosmic expansion rate. Under this extreme limit, the tightly coupled photons and baryons can be treated as a fluid in 
which the perturbations drive sound waves. The BAO effect in the cosmic microwave background (CMB) radiation 
has been subjected to extensive theoretical studies (see the early work of Peebles and Yu [42[ , a powerful analytical 
treatment by Hu and Sugiyama [43| in Fourier space, another by Bashinsky and Bertschinger [44[ in position space, 
and a review by Hu and Dodelson [IB)]). It has been confirmed and measured by CMB observations such as the 
Wilkinson Microwave Anisotropy Probe (WMAP) [46-50]. We will not discuss CMB in detail, and mainly concern 
ourselves with the aftereffect of BAO, namely its imprints on the large-scale distribution of matter. 

The imprints of BAO in the observable distribution of galaxies today was predicted in theory (see [111, [52[ , and note 
that these papers were mainly written in the language of P{k) rather than £(?")). They were first detected in SDSS 
data by Eisenstein et al. HU . In [2(} , BAO measurements were made for SDSS and 2dF survey data using the power 
spectrum, and the results were presented as a, general test of cosmological models. The usage of BAO signatures in 
the LSS as a probe of H (z) was discussed in [l(| (see also [54]-[56j|). 

The idea of using BAO scales may appear to be simple and straightforward by our description so far, but in reality 
the autocorrelation function is subjected to various distortion effects that must be accounted for. 



dP 12 = n 2 [1 + £(!•!, r 2 



)] dVidVa 



(22) 




(23) 
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First, galaxies are not comoving objects. Their apparent redshifts are inevitably a combined effect of the cosmolog- 
ical redshift and peculiar velocities (which was once contemplated by E .P. Hubble, see Section lITBl) . Peculiar motion 
distorts the apparent correlation pattern in the redshift space and makes it anisotropic (see [53, [58|). Therefore, the 
isotropic autocorrelation function £(r) fails to be a good measure. In the literature the autocorrelation is usually 
expressed as a function of scales in the radial (line-of-sight) direction tt and transverse direction a: £ = £(cr, 7r) with 
r = \J a 2 + 7r 2 . The observed £(er, tt) will be a convolution between £(r) and the peculiar velocity field. 

Second, geometry of the spacetime also distorts the correlation pattern as the observation goes into deeper distances, 
where the spacetime geometry becomes non-trivial [59]. This is not a major concern for the analyses we will review 
in the rest of this section, because the survey data were from our local section of the universe (z k, 0), and for H{z) 
measurements only some thin slices in the redshift space were used. However, future work that deals with deep survey 
data must take the geometrical distortions into analysis. 

There is also the more delicate issue of biasing, meaning that the correlation pattern of the observed "indicators" 
does not necessarily reflect that of the underlying matter distribution [60]. Among the effects contributing to the 
bias, the magnification effect by weak lensing is worthy of notice for our discussion, because it has a large effect on 
the radial autocorrelation function [6l|, [62[ . 

Using SDSS LRG samples in the redshift range 0.16 < z < 0.47, BAO signature was detected in £(a, tt) by Okumura 
et al. 63] . In their work the magnification bias by weak lensing was neglected, but in the redshift range it contributes 
little to the spherically averaged autocorrelation £o (fill ; also known as the monopole: 



1 



Zo(r) = -J a*,*)dn, (24) 

where r = V a 2 + tt 2 , and \i = it jr. In (63| the BAO peak was detected in the monopole significantly, while the 
ridge- like BAO feature was weak in the anisotropic £(<r, tt). 

Using improved LRG samples from SDSS DRs 6 and 7, and by modelling the weak lensing magnification bias, 
radial BAO detection and H{z) measurements were made in redshift slices z = 0.15 ~ 0.30 and z = 0.40 ~ 0.47 by 
Gaztanaga et al. [13] (see Figure [T] for a presentation of the BAO detection). Because these redshift slices were well 
separated, the two measurements were independent from each other. (In previous works such as [2(| the samples 
overlapped and the results at different z's were correlated.) 

These H(z) measurements were the first implementation of the radial BAO method. Due to the distortion effects, 
confirming the significance of the baryon ridge detection becomes a demanding process, since each distortion effect 
has to be carefully modelled. However, exact modelling of all the distortion effects on all scales is difficult, and when 
such modelling cannot be done exactly, these effects introduces systematic errors in the measurement of the BAO 
ridge's scale. 

Despite these, the radial BAO size method still surpasses the age method in precision. In fact, the combined 
statistical and systematic uncertainties presented an precision of ~4% in H(z) [6J, Table 3]. This is intuitively 
perceptible. As we have seen in Section UlI Al the age method is affected by the (possibly very large) systematic errors 
in age determination. Since we can measure spatial quantities of galaxies, i.e. the distribution of their positions, with 
much greater accuracy than we can do with temporal quantities related to some vaguely defined event (namely the 
time duration from star formation in the red galaxies to now), one may intuitively expect lower uncertainties from 
the radial size method than the differential age method. 

A subtle issue of possible circular logic in the analysis also contributes to the systematic errors in this method. 
In [U, a fiducial flat ACDM model and parameters were used to convert redshifts into distances, and to gauge 
the comoving BAO scales in the selected redshift slice, tbao to that of the CMB measured by 5-year WMAP, 
rwMAP = 153.3 ± 2.0Mpc (see [65]) to yield the estimation Hbao( z )'- 

Hbao(z) _ H &d (z) 



^BAO '"WMAP 



where 



H Rd ( Z ) = iW^m(l+2) 3 + (l-a„) 



and f2 m = 0.25 [6a ] . The use of a fiducial model introduces bias in all measurements, which is hard to model exactly, 
but an analysis of this effect was performed using Monte Carlo simulations so that its contribution to the systematic 
uncertainties could be assessed. The authors of [64| hence argued that the measurement results are model-independent, 
therefore is useful as a general cosmological test. The reader may also consult [20] for a different approach to this 
issue, using cubit spline fit of the distance-redshift relation so that the result could be applied to a large class of 
models without having to re-analyze the power spectra for each model to be tested. 
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Figure 1. Detection of radial (7r-direction) BAO by Gaztanaga et al. 64, Figure 13] in the full LRG sample. This is the 
correlation pattern along the 7r- direct ion, and should not be confused with the monopole pattern in Figure 3 of [6^]. The effect 
of weak lensing magnification can bee seen by comparing the solid and short dashed curves, which shows that the magnification 
systematically moves the peak location towards the higher scales. The dash-dotted (blue) curve shows the la range by allowing 
the fiducial distance-redshift relation used in the analysis to vary in a parameterized way, accounting for the systematic error 
introduced by the mere using of a fiducial model. 



A Word on the Dispute over the Radial BAO Detection. Currentl y th ere is some dispute over the claimed detection 
of radial BAO and measurement of H(z) in [64| . Miralda-Escude l67l argued against the methods in [64J and the 
statistical significance of the claimed BAO detection. Kazin et al. [68] analyzed the SDSS DR7 sample of LRGs and 
obtained similar results to [64[, but offered another interpretation using the x 2 / (degree of freedom) statistic and the 
Bayesian evidence [69j that disfavors a statistically significant detection. On the other hand, the recent research of 
Tian et al. [7(3] claims that the radial BAO feature is not a fluke, albeit certain assumptions made this re-assessment 
somewhat optimistic. The authors of (64[ also defended their work in [7l|. We refer to these variety of arguments and 
opinions to remind the reader of these ongoing investigations. Nevertheless, we believe that the general method of 
measuring H(z) using radial BAO is well-motivated and promising regardless of its current implementation, as it is 
expected to give more definitive results of radial BAO and H(z) measurement with upcoming redshift survey projects 
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Table I. The set of available observational H (z) data 



z H(z) ± la error^ References Remarks 

0.09 69 ±12 [24,^6] 

0.17 83 ±8 [26] 

0.24 79.69 ± 2.65^, [64] In the redshift slice 0.15 ~ 0.30 

0.27 77 ±14 " [26] 

0.4 95 ±17 [26] 

0.43 86.45 ± 3.68^, [64] In the redshift slice 0.40 ~ 0.47 

0.48 97 ±62 " [26] 

0.88 90 ±40 [26] 

0.9 117 ±23 [26] 

1.3 168 ±17 [26] 

1.43 177 ±18 [26] 

1.53 140 ±14 [26] 

1.75 202 ±40 [26] 



a H (z) figures are in the unit of km s 1 Mpc . 

b Including both statistical and systematic uncertainties: a = w o^ ta + cr^. 



IV. OBSERVATIONAL HUBBLE PARAMETER AS A COSMOLOGICAL TEST 



The efforts in obtaining observational H(z) data was certainly done with the goal of testing cosmological models in 
mind. In [24j the observation H(z) atzsj 0.09 was used to constrain the equation of state parameter of dark energy. 
In [28| the redshift- variability of a slow-roll scalar field dark energy potential was constrained by the differential age 
H(z) data. The same dataset was also utilized in the study of the ACDM universe, especially the summed neutrino 
masses m„, the effective number of relativistic neutrino species iV re i, the spatial curvature f2k, and the dark energy 
equation of state parameter u> [z3|. The updated H{z) data presented in [2(| was used by their authors to improve 
the results obtained in earlier papers. 

In particular, the combination of CMB and H{z) observation is a very effective way to constrain iV re i [73l see the 
reproduced Figure [5] in this paper]. In this paper we will not go further into the topic of cosmic neutrinos, which is 
intrinsically related to fundamental physics. However, we should point out a remarkable result, that the H(z) data, 
when used jointly with CMB and other late-era cosmological tests, offer valuable insight into the neutrino properties 
related to the much earlier universe, independent of Big-Bang neucleosynthesis (BBN) [74], [7f| tests. Moreover, the 
BBN constraints are obtained using Helium abundance measurements that are subjected to the systematic biasing 
effects arising from late-time neucleosysthesis. Therefore, H(z) data is an important consistency check measure in 
the presence of this systematic uncertainty [73l |. 

Figure [3] shows that adding H(z) data helps with breaking the degeneracy between spatial curvature and dark 
energy equation of state. In the ACDM universe, both the dark energy and spatial curvature becomes dominant in 
recent epochs. Therefore, separating their respective effects on the expansion of the universe becomes important, as 
well as challenging [7a, [77] . While other tests using the combination of weak lensing and B AO are likely to measure 
the curvature distinctively in the future [H|,[78j], our current knowledge of H(z) is still a valuable complement to other 
tests in the sense of DE-curvature degeneracy breaking (72| • 

The data produced by the BAO size method in 6J| is scarcer in quantity but of higher precision. In |64j they 
were extrapolated to z = to offer an independent estimation of the Hubble constant Ho, and were used to test the 
accelerated expansion of the universe. It has been demonstrated that the radial Azbao measurements is able to put 
stringent constraints over the dark energy parameters [79(. 

In the papers cited above, the parameter constraints obtained from observational H(z) data were shown to be 
consistent with other cosmological tests, such as the CMB anisotropy. In this way, the observational H(z) data 
presents themselves as a useful, independent cosmological test. In particular, it serves as a powerful tool to break the 
degeneracy between the curvature and dark energy parameters. 

These up-to-date data are summarized in Table|TJ In Figure|3]we plot the H(z) data versus the redshift. To help vi- 
sualizing the data, we also plot a spatially flat ACDM model with fi m = 0.25, Q\ = 0.75, and Ho = 72 km s -1 Mpc -1 . 

In addition to the above authors, the observational H(z) datasets have been widely used to put various cosmological 
models under test. The first adopters included Yi and Zhang [82| and Samushia and Ratra [H| who made use of the 
H(z) results of [28] in the study of dark energy. In [82j the H (z) data alone were used to constrain the parameters of 
the holographic dark energy model, especially th c parameter that determines the dynamical history of the expanding 
universe (see Figure 0. The same dataset has also been used to study modified gravity theory such as f(R) gravity 
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Figure 2. Constraint on the effective number of relativistic neutrino species, N Te \, by Stern et al. [2(| using their H(z) 
measurements by the differential age method. Dotted line plots the 5-year WMAP [8(J| likelihood, dashed line plots the 
likelihood with WMAP and Ho determined by Riess et al. [8l|. and the solid like the likelihood with WMAP, Ho and H(z) 
data. Adding H(z) data helped refining the constraint to N TC \ — 4 ± 0.5 at 1-a level. The improvement in the constraint 
by adding H(z) data is evident. Note that this figure displays the deviation of the x 2 statistic from its minimum inverted 
(Ax 2 = Xmin — X 2 )- The intersections of the Ax 2 plots with the constant A^ 2 = 4 line correspond to 2-a constraints. 



in the context of cosmology 84] . The u pda ted data in [26J and [64| have been adopted to constrain the parameters 
in more exotic dark energy models, e.g. |85l. l86j . 

Beyond parameter constraints, the observational H{z) data are also applicable in non-parametric, model-independent 
cosmological tests. For example, the Om statistic by Sahni et al. [13], defined by 

where h is the dimensionless Hubble parameter, h = H(z)/Hq. This statistic is useful as a null test of dark energy 
being a cosmological constant A, and is more robust than parameterizations of the dark energy equation of state. 
Another result for testing A that incorporates H{z) data (the £ gen test) is given by Zunckel and Clarkson 88], with the 
addition of distance information. In either paper however, the Hubble parameter data used were not the independent 
observational measurements discussed in this review, but the ones reconstructed using SNIa luminosity distances. In 
a similar fashion, it has been shown that H(z) and distance measurements can further test the spatial flatness of the 
universe, or even the Copernican Principle of large-scale homo gen eity and isotropy that is behind the mathematical 
form of the FRW metric (0) by a model-independent approach |89l. I90T] . In [§0[ the use of H(z) in some of these tests 
was demonstrated with real-world observational data reviewed here. 

Despite the wide application of the H (z) datasets in the literature, we would like to point out some issues associated 
with their usage. 

First, in some papers 85, 86, 91] that made use of H (z) data derived from radial BAO by Gaztahaga et al. {64] in x 2 
analyses, the measurement at a middle redshift z — 0.34 was used in conjunction with those from the two independent 
redshift slices near z = 0.24 and 0.43, under the tacit assumption of being independent from each other. However, this 
is not true, because the determination at the middle redshift was not made from a separate, non-overlapping redshift 



12 



0.10 
0.05 

0.00 

a -0.05 

-0.10 

-0.15 
-0.20 

-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 

w 

Figure 3. Joint constraint on the energy density corresponding to the spatial curvature, fi^, and the dark energy equation of 
state parameter, w, by Stern el al. [2(J. The large, irregular regions bounded by dark contours were from 5-year WMAP alone. 
The blue contours were obtained by adding Ho constraints. Filled regions were obtained by further adding H(z) data. The 
application of H(z) data helps with breaking the degeneracy between f2 k and w. 



slice, but from the whole sample of galaxies, including the lower and upper redshift ranges. If the data is to be used 
in quantitative works, this interdependency should not be ignored and must be explicitly analysed. A related issue is 
combining the H{z) data determined from radial BAO peaks with the Azbao data derived using the same method 
under the assumption of their independence (this practice can be found, for example, in |9l|). To be rigorous (or 
pedantic, depending on your point of view), we do not believe that this is the best way to use the data, and we insist 
on an analysis involving the (non-diagonal) covariance between these datasets. On the other hand, the combination 
of Azbao data and age-dated H(z) is mostly free from this interdependence problem, and they actually complement 
each other well [H, in particular Figures 1 and 2]. We also note that in qualitative explorations one may choose to 
relax this restriction to some reasonable extent, for example in the discussion of accelerate expansion in |64|, Section 
4.4]. 

Another topic that cold be worthy of future discussions is the possible tension between the H(z) datasets and other 
observational data. As noted by Figueroa et al. [721 ] . datasets of different physical natures and systematic effects 
can be safely combined only if they agree with each other well (see also (93[). In this regard, we note that there is 
possibly some tension between H(z) and type la supernova (SNIa) luminosity distances as shown in [92| (see Figure 
[6]). However, this apparent tension could be statistical in nature and may simply be a consequence of not having 
enough independent measurements of H{z). We hope that future expanded H(z) datasets would allow us to check 
its consistency with other data in a quantitative manner. 



V. FUTURE DIRECTIONS 

The available H(z) data have so far proven to be a useful tool in the pursuit of understanding the expansion 
history of the universe and the possible nature of dark energy. However, these datasets do not have very good redshift 
coverage. The current measurements have gone as deep as z — 1.75, and this redshift range is only sparsely covered. 
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Figure 4. Top panel — the available H(z) data from both differential age method and radial BAO size method (see Table U 
and references therein). The solid curve plots the theoretical Hubble parameter Had as a function of z from the spatially flat 
ACDM model with Q m = 0.25, CIa = 0.75, and Ho = 72 km s _1 Mpc -1 . Bottom panel — the same data, but the residuals with 
respect to the theoretical model H^d are plotted. In both panels, the z error bars on the measurements from the radial BAO 
method are used to mark the extents of the two independent redshift slices in which the BAO peaks were measured. 
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Figure 5. Parameter constraints for the holographic dark energy model in the f2 m -c plane, by Yi and Zhang [82] . The 
constraints were obtained using age-determined H (z) data in [2|| alone. The cross in the lower-left marks the best-fit value. 
The dash-dotted, solid, and dotted contours marks the 68.3%, 95.4%, and 99.7% confidence regions respectively. Although 
some degeneracy exists, it is evident that the data favor models with c < 1. 

There is also another issue of the large error bars associated with the H(z) figures from the differential age method. 
On the other hand, the collection of more and higher quality H(z) data will not only help us constrain the parameters, 
but will also allow us to understand the possible tension between H(z) and other cosmological tests. The latter is 
important, because tension is usually an indicator of systematic errors in the data. By understanding the tension, we 
may finally conquer the systematic effects that have not yet been modelled well enough. 

In this section, we will describe a few directions of future cosmological observations and their implications in the 
measurements of the Hubble parameter. 



A. Future Improvements in the Differential Age Method 

The relatively large uncertainties in the differential age method could be partially compensated if future datasets 
could offer better coverage in the redshift range accessible by this method. Using mock data, we recently estimated 
that future H(z) datasets would offer similar or even higher parameter-constraining power compared with current 
SNIa datasets if it could add as many as ~60 independent measurements to cover the redshift range < z < 2 [95| . 
To achieve this level of data coverage, future surveys must be able to offer a large sample of LRGs to be used in 
age-dating. According to [28[ , the Atacama Cosmology Telescope (ACT) (9(| can be utilized in the future to identify 
passively evolving, red galaxies by their Sunyaev-Zel'dovich effect. These galaxies can in turn be spectroscopically 
measured and age-dated, and it has been estimated that they could yield ^1000 H(z) measurements. This means the 
quality of current differential age H (z) data can be expected to increase sig nificantly. 

The error model used in the analysis of differential age H(z) data in [95| was empirical, which may have underes- 
timated possible future improvements. In [36| it has been estimated that H (z) may be measured within 3% relative 
error at z 0.42 in realistic observations if the star formation systematics could be properly accounted for. This level 
of precision is on par with the current status of the radial BAO method, and we hope it could be achieved in the near 
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Figure 6. Possible tension between H(z) and type la supernovae data depicted in the \ 2 fitting results for the spatially flat 
XCDM model (similar to ACDM, except that the dark energy equation of state parameter cj is set free instead of being fixed 
at to — —1). The SN data favor a phantom dark energy with to < —1 while other data, including observational H(z) (OHD), 
are consistent with ACDM. The OHD used in this figure were the measurements by [28j] using the differential age method, 
and the SN data were from [93]. The RBAO contours were found using the Azbao data in [63]. Confidence regions are 68.3%, 
95.4%, and 99.7% respectively. This figure first appeared in [92l . Fig. 4]. 



future. 



B. Future Improvements in the Radial BAO Size Method 



The radial BAO size method has already been demonstrated to provide highly accurate H(z) measurements. How- 
ever, this accuracy came at a cost, for spectroscopic data must be taken for the great number of galaxies under survey 
to find their redshifts, which is time-consuming. Fortunately it turns out that for low redshift ranges, photometric 
redshift surveys can be a sufficient and promising approach [97H99I ] to the detection and measurements of radial BAO 
features in the autocorrelation function. Photometry has several advantages over spectroscopy - it is cheaper, faster, 
and able to reach fainter sources. 

Shortly be fore this review is written, the WiggleZ redshift survey [lOOj j of emission-line galaxies produced its first 
data release [lOl| . As the data is being released, it is expected that the radial BAO signal can be put to further 
scrutiny (68j . 

The BAO method is unique in that it allows us to reconstruct the cosmic expansion through a vast range of eras. 
Unlike the differential age method in which the observable indicators of time are located within a limited redshift 
range, BAO signal detection is possible as long as the distribution of matter, regardless of its form, can be traced. 
Even if the current implementation of the radial BAO metho d is mainly confined in the redshift range of z ps 0, future 
redshift surveys such as the planned SDSS III project jlQ2| ] are designed to reach into deeper universe and me asure 
H(z) at redshifts up to z ~ 2.5 by observing the Lyman-a forest absorption spectra of high-redshift quasars (see pi 
for a discussion of high-z measurement of radial BAO and H(z) and its implication for dark energ y, an d pill05j for 



i nergy, an c 

numerical simulation studies). Recently, in the wake of the proposed Euclid satellite project |l06l Il07| . the enormous 
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potential of space-based redshift surveys in the determination of H{z) and other parameters has been studied in [108| . 
Finally, the proposed observational programs of the 21 cm backgrou nd may f urther extend our knowledge of H(z) 
into even deeper redshift ranges before or near the reionization era the "dark ages" that have not been 

extensively explored by current observations yet. 

It is also worth noting that the previous works on the analysis and measurement of H(z) from the clustering of 
LSS have mostly concentrated on the BAO features alone. However, Shoji et al. [l!2j | shows that accurate estimates 
of H(z) and Da{z) could be made using the full galaxy power spectrum in the extraction of cosmological information 
instead of BAO features alone, provided that the non-linear clustering effects are well controlled. We hope that the 
future redshift surveys observations, as well as advances in better understanding of nonlinear-regime redshift-space 
distortions, could lead to successful realization of their method. 



VI. SUMMARY 



In this paper, we reviewed the current status of observationally measured Hubble parameter data. We presented 
the principle ideas behind the two important and independent methods of H(z) measurement, namely the differential 
age method and the radial BAO size method. Both methods have been successfully implemented over the years to 
yield H(z) data that are of varying precision and redshift coverage, and the up-to-date results have been summarized 
in Table HI These data are valuable for the study of the expanding universe. They have seen wide application by 
cosmologists to put various cosmological models under test, and to constrain important cosmological parameters either 
independently or in conjunction with data of different physical natures. However, we also pointed out several issues 
in the usage of observational H(z) data. Finally, despite some current shortcomings, we find the H(z) data of great 
potential, as future observational programs can be expected to improve significantly the quality of H(z) data that 
may lead us into unexplored realms of the universe. 
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